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Abstract 

A smooth unclustered dark matter component with negative presure could rec- 
oncile a flat universe with the many observations that find a density in ordinary, 
clustered matter well below the critical density and also explain the recent high red- 
shift supernova data suggesting that the expansion of the universe is now accelerating. 
For a perfect fluid negative presure leads to instabilities that are most severe on the 
shortest scales. However, if instead the dark matter is a solid, with an elastic resis- 
tance to pure shear deformations, an equation of state with negative presure can avoid 
these short wavelength instabilities. Such a solid may arise as the result of different 
kinds of microphysics. Two possible candidates for a solid dark matter component are 
a frustrated network of non-Abelian cosmic strings or a frustrated network of domain 
walls. If these networks settle down to an equilibrium configuration that gets carried 
along and stretched by the Hubble flow, equations of state result with w = —1/3 
and w = —2/3, respectively. One expects the sound speeds for the solid dark mat- 
ter component to comprise an appreciable fraction of the speed of light. Therefore, 
the solid dark matter does not cluster, except on the very largest scales, accessible 
only through observing the large-angle CMB anisotropy. In this paper we develop a 
generally-covariant, continuum description for the dynamics of a solid dark matter 
component. We derive the evolution equations for the cosmological perturbations 
in a flat universe with CDM+ (solid) and compute the resulting large- angle CMB 
anisotropy. The formalism presented here applies to any generalized dark matter 
with negative pressure and a non-dissipative resistance to shear. 
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1. Introduction 



Most determinations of the cosmological density parameter Qq = (p/p C rit), where 
Pcrit = (3/8irG) Hq, now indicate that Q m ~ 0.2 ± 0.1, a value well below the 
Q m = 1 value suggested by flat cosmological models. (For a nice review of the current 
observations see ref. 1.) Most of these techniques for determining Vt m , however, are 
sensitive only to matter that is clustered gravitationally and do not rule out a smooth, 
unclustered component that could make up the difference between the observed value 
of fl m and unity. 

The earliest proposal for a smooth component is the cosmological constant A, first 
introduced by Einstein, 141 later denounced by him,' 51 and more recently resurrected 
to reconcile the observations with a flat universe.' 61 The cosmological constant is 
somewhat of an embarrassment for theoretical physics because dimensional arguments 
would suggest A « M p ; 4 , a value more than a hundred orders of magnitude too 
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big! Perhaps some not yet discovered symmetry makes A vanish exactly, but at 
this point we lack even the vaguest idea of what kind of symmetry could do the 
job. Supersymmetry somewhat mitigates the difficulty, making M$sb 4 rather than 
M p i 4 the naive guess for A, but even with supersymmetry if A does not vanish, a 
formidable fine tuning problem persists. A large £7a overpredicts the number of 
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gravitationally lensed quasars. As an alternative to A, it has been proposed that 
there could exist a very light, extremely weakly coupled scalar field that could act as 
a temporary cosmological constant, even though the true value of the cosmological 
constant vanishes exactly' 251 But this requires a particle of implausibly small mass, 
somewhere in the neighborhood of 10 -33 eV. 

In this paper we discuss another possibility: a solid dark matter component with 
significant negative pressure. Here significant means that the negative pressure, or 
equivalently tension, of the solid matter component is comparable in magnitude to its 
energy density. An equation of state with large negative pressures can lead to sound 
speeds comparable to the speed of light, so that the Jeans length for this component 
is enormous, comparable to the size of our present horizon. Consequently, the solid 
dark matter component does not cluster except on extremely large scales. Because 
of this the low measurements of Q can be reconciled with a spatially flat universe. 
The clustering of the solid dark matter component on very large scales is accessible 
to observation only through its effect on the large-angle CMB anisotropy. 

A solid dark matter component can also help explain the the recent observations 
of distant Type la supernovae that suggest that the universe is now accelerating.' 9,101 
For the expansion of the universe to accelerate some exotic form of matter with 
w = (p/p) < —1/3 is required. A perfect fluid with w < is not possible because its 



1 



sound speed would be imaginary, indicating instabilities whose growth rate diverges 
as the wavelength approaches zero. Such a fluid would clearly be unphysical. For a 
solid, however, real sound speeds are possible because a shear modulus of sufficient 
magnitude removes these instabilities. 

In this paper we explore the dynamics of a solid dark matter component by devel- 
oping a continuum description for such a component within the framework of general 
relativity and incorporating the solid dark matter component into the linearized the- 
ory for the evolution of cosmological perturbations. In particular we explore the 
consequences of such a component for the large angle CMB anisotropy. 

A solid dark matter component could arise from a variety of different micro- 
physics. Two known ways such a component could arise are from networks of frus- 
trated cosmic strings 119 ' 18 ' 171 or domain walls! 24 ' 231 While the simplest Abelian cosmic 
strings obey a scaling solution so that the number of strings per horizon volume re- 
mains constant, for non-Abelian cosmic strings topological obstructions prevent the 
intercommuting necessary for the breakup long strings that leads to scaling behavior. 
The nonunit elements of the fundamental group 7r\(G/H) classify the possible types 
of cosmic strings. When two strings whose windings or magnetic fluxes are described 
by non-commuting elements of tt\{G/H) try to cross, the strings cannot pass through 
each other without forming a third string between them. This has the effect of pre- 
venting crossings because the tension of the intermediate strings tries to pull the two 
strings back toward their previous uncrossed positions. It is possible that these effects 
lead to a scaling solution albeit one with many more strings per horizon volume, but 
the simulations by Pen and Spergel suggest that the strings settle down to a stable 
configuration which subsequently gets carried along with the Hubble flow. In a forth- 
coming paper, we show that stable string configurations do exist which strengthens 
the case for a string dominated universe. Similar simulations of domain walls by 
Kubotani suggest domain walls in models with many types of domain walls exhibit 
similar behavior. A cellular foam type structure in equilibrium forms with several wall 
meeting at linelike junctions. A string-dominated universe gives w = —1/3, which in 
the absence of other dark matter would give a universe that is neither accelerating 
or decelerating. A symmetry breaking scale of a few TeV and a string separation 
today of a few A.U. would give ^string today in the interesting range between zero 
and one. (This estimate is subject to considerable uncertainty because a range of 
string density at formation is possible and the length of the transient before string 
dominated behavior takes over is uncertain and model dependent). For domain walls 
carried along with the Hubble flow, w = —2/3 and a symmetry breaking scale of a 
few MeV and a mean domain wall separation of some tens of parsecs are suggested 
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(subject to the same uncertainties). The fact that new physics occurs at larger energy 
scales than for quintessence or A is a positive feature of these scenarios. It should be 
stressed that the formalism in the paper applies equally well to a solid component 
with the same continuum description but of completely different origin. 

A slightly different type of solid dark matter has been proposed by Eichler to 
explain certain aspects of large-scale structure.' 321 In this scenario a solid condenses 
and subsequently fractures when stretched beyond its breaking point by expansion 
of the universe. The dark matter contemplated in this paper does not fracture. It 
can experience unlimited stretching without becoming in any sense weakened. For 
a solid composed of frustrated topological defects it is easy to see why ruptures or 
fractures do not occur. The constituent defects lack a preferred size. Upon stretching 
or shrinking, the transverse structure of the defects remains unchanged. This is quite 
unlike an ordinary solid composed of atoms, for which quantum mechanics establishes 
a preferred length for the chemical bonds. 

The organization of this paper is as follows. In Sect. 2 of this paper we develop a 
generally covariant description of the dynamics of a continuous medium (such as the 
string network) in curved space. For the spacetime with the metric G^ v = a 2 (i]) + 
h^] where is regarded as small, we expand the action to quadratic order and 
compute the equations of motion and the resulting stress-energy for the solid dark 
matter component. In Sect. 3 we combine the results of Sect. 2 with the linearized 
theory of cosmological perturbations using Newtonian gauge and derive the equations 
of motion for a spatially flat universe with cold dark matter (CDM) and a solid dark 
matter component. In Sect. 4 compute the large- angle CMB anisotropy for models 
with SDM. Finally, in Sect. 5 we present some concluding remarks. In this paper we 
use the sign convention rj^ = diag[— 1, +1, +1, +1]. 

2. Continuum Description 

This section presents an action that describes the dynamics of a dissipationless 
elastic medium in curved space. Although developed to describe the response of a non- 
Abelian string network to metric cosmological perturbations, this formalism applies 
to other forms of solid dark matter and to a wide variety of situations involving 
continuous media in curved space. The problem of describing the dynamics of a solid 
within the framework of general relativity has been previously considered by Carter 
and Quintana in the study of the crusts of neutron stars and others. In this 
section we present a self contained treatment particularly suited to the consideration 
of linearized perturbations in an expanding universe. 

A continuous medium is a kind of three-dimensional membrane but quite dif- 
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ferent from the now much discussed 'fundamental' branes. A continuous medium 
has internal structure. As the medium moves, each constituent particle traces out 
its own worldline. For continuous media the allowed reparameterizations are limited 
to reparameterizations that preserve these worldlines. Geometrically a continuous 
medium may be regarded as a congruence of worldlines. We use a three dimensional 
coordinate y to parameterize these worldlines. The coordinate t is an arbitrary time 
coordinate parameterizing the direction along these worldlines. 

A metric h ab is defined on y-space. The volume induced by this metric indi- 
cates the density of worldlines and the additional internal metric structure provides 
a reference with respect to which to characterize pure shear (i.e., volume preserving) 
deformations. In the classical exposition of elasticity theory (e.g. as in the book by 
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Landau and Lifshitz ) h ab = 5 ab = (constant) . However, when a solid is formed in 
a warped (i.e., curved) spacetime, there generally does not exist any coordinatization 
of the worldlines so that h ab = 8 ab . For the moment let us imagine ourselves in the 
instantaneous rest frame of a volume element of the medium, choosing t so that d/dt 
is orthogonal to d/dy a (a = 1,2,3). The background spacetime metric G^ v induces 
the following metric on y-space: 



dX^dX^ 

9ab - dya dyb . (2.1) 



For an arbitrary time parameterization, where d/dt is not necessarily orthogonal to 
d/dy a , the induced metric may rewritten as 

9ab - G, v dya dyb (2.2) 



where G$ = U^U V + G^ and U» = ^/y^^T^T- G $ projects out displace- 
ments along d/dt. 

The local deformation state of the medium is determined by comparing g ab to 
h-ab — by the three principal values Ai,A2,A3 of g ab with respect to h ab . (We assume 
that the medium is isotropic, for otherwise more structure than h ab alone is required 
to characterize the deformation state of the medium.) In ordinary elasticity theory, 
h a b — &abi 9ab 1S the strain tensor, and Ai, A2, A3 are its principal values. The 
scalar invariants = h ab g ab , # (2) = h ab g bc h cd g da , # (3) = h ab g bc h cd g de h e f g fa suffice 
to characterize completely the principal values, and the local density in the local 
instantaneous rest frame with respect to the h volume element may be expressed 
as P(h) = P(h)(9(i),9(2),9(3))- In terms of the principal values g^ = Ai + A 2 + A3, 
9(2) = Ai 2 + A 2 2 + A3 2 , and g {3) = Ai 3 + A 2 3 + A 3 3 . 



4 



It follows that the action is 



/r J dX^ dX v 

dt J d 3 y Vh p(g {1) ,g {2 ), 9(3)) y ~ G ^-Q t — gf- ( 2 - 3 ) 

One may view (2.3) as a generalization of the free particle action S = —m J dr. 
If p(g(i), <7(2), ^(3)) = (constant), the action (2.3) merely describes a congruence of 
non-interacting, freely-falling particles. However, in the general case the co-moving 
density with respect to the internal coordinates varies as deformations of the medium 
alter its internal energy. The potential energy of the medium resides in the function 

P(9{i),9{2),9{3))- 

We now recast the action (2.3) into a more familiar form by considering an elastic 
medium in almost flat space. We assume a spacetime metric G^ v = rj^, + h^, an 
internal metric h a \y = 5 a b + b a b, and X % = y % + £*(y, i) where h^, b a i,, and C{y,t) are 
regarded as small. We also set X° = t. 

We expand 

P = Ps + r s (^j + K s 2 + p. S {5)ab S$ y (2.4) 

Here r s is the tension and K s and /i s are the compressional and shear moduli, respec- 
tively, and the tensor is the pure shear component of the strain tensor. From 
the relation (valid to quadratic order) 



$V\ _ y/\(jab\ _ y/\5ab + S a b\ _ 1 + § + %" 



2 o . . M 

a 



it follows that 



6V s s 2 SijsV b h.W b 2 bs 

o + ^ + ^-^- ( 2 -6) 



V 2 8 4 2 4 8 1 
Here to quadratic order 



Safe = {S\ + C , a (y, t)}{*y + hy(X) + (& + h 0l )fe + h ,)}{<F 6 + ^(y, t) } - 5 a fe 
= £a,b + £b,a + h a fe + ^Vfc h a& + &, a &,& + h a j£ i;6 + hfo^a + £ a 4 
+ hoa^b + h b^a + hOahofe, 

(2.7) 
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and to quadratic order 



£,a,a + n aa + V c h aa £,a,b£,b,a £a,a£&,6 'Ca,a^&& Macfcbb 



V 



+ gh aa h 6& - -h afe h a& + -£ a ,ahbb - -jbaa + -^hbhb + gbaahb + T^a ( 2 ' 8 ) 
+ hoa'Ca + 2 h 0ah()6- 



To linear order (which is sufficient here), the pure shear part of the strain is 



&ab = £a,b + ib,a + Kb ~ b a b ~ ^ab ( 2 £ c , c + h cc - b cc ) . 



1 



(2.9) 



Combining and expanding to quadratic order, we obtain 



dt / d 6 y 



, 6 6 2 6 a6 6 a6n 

H 1 — 

2 8 4 



x 



1, , 1 1 1. 1, 1, , 1, , 

- -^haabbb + gh aa h bb - -h ab n a& + -£a,ah&6 - -o aa + -o ab b ab + -o aa o b b 



+ 



+ hoa'Ca + ^hoahofr j 

r 1 • 1 • "i 2 

+ Psj^a.fo + £&,a + ha& - &afc - ^ab (2 £c,c + h cc - b cc ) j 



(2.10) 



X 



1 • 



z 2 



and expanding out to quadratic order and omitting terms that do not contribute to 



6 



the equations of motion or to the stress-energy T^, we obtain 



= /*/ 



Ps - T s 



s \ <-2 



C + ( Ps - T S ) hQiC 



r noo 



h o b aa h Q h, 



00 



-^e c V c )h 00 } 



1. 1,1 1 1 1-1 

+ 2^,0 hfefc + -h o o aa - -h o h aa + -h aa h bb - -h ab h ab + -h 0a hoa J 

-^{r„ + ^.-^,} 2 

- /^s |£a,6 + £b,a + ^afc ~ b ab - -5 a& ( 2 £ C)C + h cc - b cc ) j 



(2.11) 



It follows that the equation of motion for £ is 



(p s - T s )^i + (p s - r s )hoi - -(ps — r s )Vihoo - K s 2£jji + hjjj - bjj^ 



- Ps 



444 



= 0. 



(2.12) 



With the coupling to gravity ignored, which is a reasonable approximation for short 
wavelengths, the sound speeds 



9 o ps 2K S 

c S = 



(Ps - Ts) 



Cy 



^Ps 



(Ps ~ T s ) 



(2.13) 



follow for the longitudinal (scalar) mode and the two transverse (vector) modes, 
respectively. 

We next compute the stress-energy, defined by the relation 



5s= \S d " x ^ t>iu (5G ^ } 



(2.14) 



We use the expansion l/y/-\Gp V \ = [l + hoo/2 - ha/2 + 0(h 2 )] and the relation 
A w 2 aa + B(w ab - \5 ab Wcc) 2 = (A - \B) w aa 2 + B w ab w ab . It follows that (to linear 
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order) 



T ° = 


~Ps 


- \{ps - r s ){bu 


A i — 


(Ps ' 


-T S ) (C + hoi), 


1 I — 


-T a 


Sij 




- 2 


K s - ^/i s J 6^ 



(2.15) 

We now turn to the continuum description of the solid dark matter component 
with the expansion of the universe included using conformal time 77 so that the back- 
ground metric becomes 

G^ = a 2 { V )- [77^ + V]. (2.16) 
We modify the notation for the expansion of the co-moving density as follows 

P = Psa\r ] )ll + f s (^j +K s (y) + p s S {5)ab S (5) A , (2.17) 



so that f s , K s , and p, s are dimensionless. We take the exponent A to be constant, 
although the generalization of this is straightforward. The physical density scales as 
Pphys — Pcm/a?(v) = Ps a x ^(r/). Given A, the dimensionless parameters 



f. = * k, = t^JL (2.18) 
3 18 V ; 



are fixed, and j2 s is variable only within the range 



max 



0,— K s 

' 8 



<p s <^—^, (2.19) 



obtained by requiring stability on short wavelengths and subluminal longitudinal and 
transverse sound speeds. For frustrated non-Abelian strings, f s = 1/3 and K s = 
— 1/9; for a frustrated network of domain walls, f s = 2/3 and K s = —1/9. 
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For the expanding universe, the action (2.10) is modified as follows: 



S=Jdr ] Jd 3 y p s a x+1 (r)) 



b b 2 b ab b abl 

H 1 — 

2 8 4 



I. 1 1 1 1 1 1 

- -^h-aabbb + gh aa h &fe - -h afe h a& + -4 a , a h 6 5 - -0 aa + -b ab b ab + -b aa b bb 
f 1 V 



(2.20) 



X 



The equations of motion are modified to 



(1 - fs) [ii + (A + + (1 - f 8 ) [ha + (A + l)-hoi] - -(1 - f s )Vih 00 

4. .. 4, 4 



2^+h 



/is 



In the expanding universe, the stress-energy (to linear order) is 



= 0. 
(2.21) 



T ° = a^p s 
Tj> = al x -Qp, (l-r s ) (e + h Qi ) 



1 - -(1 - T s )(bu - ha - 2£ M ) 



-t» Aij - '2 [K a - j 5 ,j Uo,o + ^h„„ - hj aa 



(2.22) 
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3. Coupling to Gravity 

We choose Newtonian gauge (which is equivalent to the gauge invariant formalism 
of Bardeen), so that the metric takes the form 

ds 2 = a 2 (r]) \-drj 2 (l + 20) + dx i dx j {<% (l - 2^) + hg° + } ] (3.1) 

(V) (T) 

where • and • are pure vector and tensor parts of the spatial-spatial metric 
perturbation, respectively. In dealing with cosmological perturbations it is convenient 
to define any vector or tensor that can be expressed by taking derivatives of a scalar 
quantity as scalar. Likewise, a tensor that can be expressed as a derivative of a vector 
is regarded as vector. With these definitions the linearized equations separate into 
independent scalar, vector, and tensor blocks. We assume a flat universe filled with 
cold dark matter and a solid dark matter component. 
We decompose the displacement field 

t = tW+t<Y) t (3.2) 
and the internal metric of the solid dark matter component 

hj = 26<? i v + 6<S (M^iM!) + 6 (V> + 6 (J). ,3.3) 



3.1. Scalar Perturbations 

We first consider the scalar perturbations. The linearized Einstein equations for 
the scalar sector are: 



vV - m{i) + n<j>) 

8ttG) [ -p c 5 + O ° ] , 

1 1 r • Ho 
a* l 

(8ttG) 



(3.4a) 



PcV\ 



(3.46) 



5G (S-tr)J = 4 



i> + 2Hijj + H<j) + {2H + Ti 2 )<P + -V 2 (0 - V) 



Si 3 
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(8ttG9 e^-^V, 



5G (S-ntr)3 



(3.4c) 



(3.40 



where tr and ntr denote the pure trace and traceless scalar parts of the spatial-spatial 
tensors, respectively. The dots represent derivatives with respect to conformal time, 
Tt = (a/a), p c = Q c dm (3/87rG) H 2 a~ 2 , and V{ = vu (i.e., v is the potential for the 
scalar part of the velocity field). The covariantly divergenceless tensor is the 
perturbation in the stress-energy of the solid dark matter component. 
Eqns. (3.4a) and (3.4b) may be combined to obtain 



2 



^cdm ( bcdm - 3H v cdm \ H f-@o° - 37iv) 

> Pcrit V 



(3.5) 



where 9° 'i(string-scalar) = V|j (i.e., Vis the potential used to represent the scalar com- 
ponent of the solid dark matter component momentum density). Similarly, eqn. (3.4d) 
may be recast as 



v*v J ' - -5i j v 2 (v> - 0) = m 



Q(S-ntr) P 



Pcrit 

The equations of motion for the cold dark matter (CDM) component are 

<W = -(V • v cdm ) + 3ip, 



(3.6) 



(3.7) 



For the scalar mode of the solid dark matter component, the equation of motion is 



(1 - r 8 ){tW + (A + l)H + V0} = (<2k s + 
Finally, we have the equations 

_ -{Ps-Ts) 



16 p s 



(3.8) 







Q(S-ntr).J = _ 4 ~ 



a( 3 - A ) 



1 1 

-^ + 3^ + 3^ 



(3.9) 



b (S-ntr) 



S a(3-A) 

With the equations of motion for the cosmological perturbations including the 
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solid dark matter component derived, we now turn to initial conditions. For each 
wavenumber k there exist four modes: two growing modes and two decaying modes. 
This pair of perturbations corresponds to the two distinct ways in which the solid 
dark matter generates and alters the growth of perturbations. If prior to the phase 
transition that produced the solid dark matter, there were pre-existing curvature 
fluctuations, then the presence of the solid dark matter alters the evolution of the 
fluctuations. This is similar to the way in which a cosmological constant, neutrinos, 
or quintessence alters the evolution of the fluctuations. In addition, the generation of 
the solid dark matter (e.g., string formation) at a phase transition can produce new 
fluctuations. These entropy fluctuations correspond to variations in the dark matter 
density at the surface T = T pt , where T pt is the phase transition temperature. These 
white noise entropy fluctuations are likely to be small on scales large compared to the 
horizon scale at the phase transition. 

We focus on the effect of the solid dark matter on the evolution of pre-existing 
scalar, vector and tensor fluctuations. Inflation generates primarily scalar and tensor 
fluctuations; however, we include the vector term for completeness. On the surface 
T = Tpt the solid dark matter component inherits as its intrinsic spatial metric the 
metric on this surface induced by the background spacetime metric. Specifically, for 
small perturbations this gives 



Kb 



-2ij + 2 ()n " 1 



3 (1 + w rad ) 



Kb, (3.10) 



or, equivalently, 



b (S) - f I 1 6r b {S) -0 (3 11) 

K - y j + 3{l + Wr y b ntr-U- (3.11) 



Initially, we assume that £ = £ = 0. The second term in (3.10) arises from the shift in 
time of the surface of constant density relative to the surfaces of constant cosmic time 
for Newtonian gauge. Since the wavelengths of interest at this point lie far beyond the 
Hubble length, we have ignored perturbations in the velocity of the radiation fluid. 
We assume that (krj) <C 1 and that only the growing mode is relevant. 

The perfect fluid analogue of the above is as follows. For temperatures T > T pt , 
the universe is filled with a single perfect fluid, which at T = T pt instantaneously splits 
into several uncoupled perfect fluid components, labeled by (A — 1, . . . , N). In this 
case the matching condition is 5/(1 + w) = 5 a/ '(1 + wa) for all A and all velocities 
may be neglected. While the Lagrangian formalism developed in this paper rather 
than the more familiar Eulerian formalism could be used to describe this perfect fluid 
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situation, for non-Abelian strings, and similarly any solid with harmonic resistance 
to shear, the more general Lagrangian description is required. 

We now consider the subsequent evolution given these initial conditions. In the 
situations of interest the solid dark matter component is formed well before radiation- 
matter equality and the solid dark matter component contributes negligibly to the 
matter density of the universe compared to other components until well into the 
matter dominated epoch. 

Assuming either complete matter or complete radiation domination gives 5 = 

( S) 

—2(f) = —2ip on superhorizon scales. Consequently, b tr = — (3/2)ip during radiation 
domination on the scales of interest. During the matter-radiation transition ip drops 
by a factor of (9/10) while b tr does not change; therefore, during matter domination 
on superhorizon scales b\f^ = —(5/3)ip. 

To follow the evolution of the perturbations through the transition from matter 
to solid dark matter domination, it is convenient to define the variable S — (V • £). 
The equation of motion becomes 




(3.12) 



and the sources become 




(3.13) 



It follows that 




(3.14) 



and 
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The evolution of the gravitational potentials after the radiation epoch may be 
computed by solving equations (3.12), (3.14), and (3.15) combined with the initial 
conditions S — S — 0, tp — ipinit, bffi = —(5/3) ipinit, and ip — 0. For a non-Abelian 
string network stretched by the Hubble flow, the evolution of the scale factor is given 
by a(rj) = a [cosh?] — 1] where fl m = sech 2 [?7/2], just as for a hyperbolic universe with 
CDM. 

As a practical matter, it is better to use synchronous gauge to evolve the per- 
turbations because synchronous gauge is better behaved on superhorizon scales. In 
Newtonian gauge the density and velocity perturbations on superhorizon scales are 
not small. This results from the warping of surfaces of constant cosmic time required 
to make the spatial part of the metric conformally flat. As a consequence using the 
constraint equations to determine the Newtonian potentials and their time derivatives 
involves delicate cancellations between large quantities, cancellations that become in- 
creasingly delicate as one passes to earlier times. Synchronous gauge, on the other 
hand, is more robust. For the adiabatic growing mode synchronous gauge (with the 
amplitude of the gauge modes set to zero) rapidly approaches a co-moving, constant 
density gauge as one passes to superhorizon scales. With the convention' 261 

ds 2 = -di] 2 + a 2 {rf) 8^ + h(k, rj)kikj + 6T/(k, rf) {kjtj - ^ij^j dx l dx j , (3.16) 

at early times during the radiation epoch on superhorizon scales one has 5 ~ O(fcV) 
and 9 ~ 0(/c 4 r/ 3 ) for all components contributing to the stress-energy and h ~ 
0(k 2 r] 2 ), and the only appreciable perturbation is rj ~ 0(1). For the solid dark 
matter for initial conditions on superhorizon scales it is an adequate approximation 
to set £ = ^ = and bij = 6rf(k,r]) — The Newtonian potentials may be 

calculated from the synchronous variables according to 



(h + 6ri)+H(h + 6r])\, 



(3.17) 
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3.2. Vector Perturbations 



For completeness in this subsection we give the evolution equations for the vector 
sector. Although for each wavenumber k the solid dark matter component has two 
dynamical vector degrees of freedom, for the usual inflationary models combined 
with a solid dark matter component these vector modes are not excited. As before, 
on the initial surface at T = T p t the solid dark matter component inherits as its 
intrinsic metric the metric on this surface induced by the spacetime metric, but if the 
hjj^ = 0, it follows that = 0. Similarly, on this surface on superhorizon scales 

£( v ) = £( v ) = 0. We also have the gauge condition = 0. 
For the CDM the equation of motion for the vector modes is 



cdm ' rL cdm 



-V-h 



(V) 



(3.18) 



Similarly, for the two vector modes of the solid dark matter component, the equations 
of motion are 

(1 _ fs) ^(V) + (1 + x)n £V)} = + V • - V • b^}. (3.19) 

The Einstein equations for the vector sector are 



- V,- • iSP 



t 2 v 3 



a 



(V) r^(V) 

KG) pcdmVcdm i + ©Qi 



(3.20) 



Finally, we have the equations 



en 



.(3-A) 



:i-r.)6°°, 



900/ = 



(3.21) 
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3.3. Tensor Perturbations 



Primordial tensor perturbations, such as those generated during inflation, are 
influenced by the solid dark matter component. 

The linearized Einstein equation for the tensor sector is 



1 



(T) 



72i;(n 



'*:) 



V 
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The tensor stress-energy from the solid dark matter component is given by 



e 



— 4/X s Pcrit 



which may be inserted into (3.22) to obtain 

hjf + m h { P - v 2 hJP + i2^ sdm 7Y 2 



h«r-h«j» h =o) 



0. 



(3.22) 



(3.23) 



(3.24) 



Physically, the response of the solid dark matter component contributes a mass term 
to the gravity waves. 

4. Implications for the CMB Anisotropy 

We now explore the consequences of a solid dark matter component for the pre- 
dicted CMB anisotropy. The scalar contribution is given by the Sachs- Wolfe formula 



AT 



Vo 

(M) = ^(n„M,i».) + Jdr, (^ + |^ 



Vis 

Expanding the CMB anisotropy into spherical harmonics 

AT 

lm 



(4.1) 



r=Vo~V- 



T 



-(M)=X>imlWM), 



(4.2) 



one obtains the following expression for the expected multipole moment variance: 



Vo 



-1 2 



-<S>(r}ls;k)ji(kri s )+ / drj ji(kr) { ^(77, k) + $(77, k) 



Vis 



ci = (\a lm \ z ) = J dkP(k) 


(4.3) 

Here P(k) is the primordial power spectrum, which we set to P(k) oc 1/k, indicat- 
ing a featureless, scale-invariant, Harrison-Zeldovich primordial spectrum. Here the 



16 



functions ^(rj, k) and $(77, k) indicate the time dependence of the growing modes of 
wavenumber k of ip and 0, respectively, and are normalized to unity as rj — > 0. 

On small angular scales (£ > 30), the solid dark matter component does not play 
a significant role in determining the CMB anisotropy because on these scales the 
anisotropy is almost exclusively determined by what happens at the surface of last 
scatter when Q so iid is negligible, and the contribution of the integrated Sachs- Wolfe 
term on these scales is negligible. On larger angular scales, however, the contribution 
through the integrated Sachs- Wolfe term of the decay of the gravitational potential 
contributes significantly to the low-£ moments, and since the details of how the po- 
tential decays depend on the dynamics of the smooth component, one expects the 
dynamics of the solid dark matter component to play an important role in determining 
these CMB moments. 

To illustrate the effect of a solid dark matter component, we compare the evolu- 
tion of the gravitational potentials and the large-angle CMB moments for the follow- 
ing nine cosmological models, some with a solid dark matter component and others 
included for purposes of comparison: 

0. A flat Q m = 1 universe. In this critical universe filled with cold dark 
matter the gravitational potentials and ip remain constant at late times, so 
there is no integrated Sachs- Wolfe contribution. 

1. A hyperbolic VL m = 0-3 universe. The scale factor a (rj) for this hyper- 
bolic universe with no dark matter other than a subcritical cold dark matter 
component evolves exactly as the scale factor for the string dominated uni- 
verse. However, the negative spatial curvature and differing decay of the 
gravitational potential at late times leads to a different shape for the low-£ 
CMB moments. 

2. A flat Q m = 0.3 string dominated universe. In this flat universe with 
a subcritical density of cold dark matter and the remainder in a network of 
frustrated non-Abelian strings, the physical density in the solid string network 
component scales as p s oc a -2 , becoming the dominant form of matter at late 
times. The behavior of the solid string component depends on the strength 
of the resistance to pure shear jl s . We consider the following three subcases: 

2a. cs = 0, cy = 1/2. (fl s = 1/24). 

2b. c,<? = I/a/3, cy = 1/V2. (fi a = 1/12). 

2c. cs = 1, cy = 1. (p, s = 1/6). 

3. A flat Q rn = 0.3 domain wall dominated universe. In this flat cosmo- 
logical model a network of frustrated domain walls formed in a late-time 
phase transition gives a density that scales as p s = 1/a. As for the string- 
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dominated universe, we consider three subcases: 

3a. C£ = 0, c v = 1/V2. (fl s = 1/24). 

3b. cs = 1/V3, cy = a/3/2, ijxs = 3/48). 

3c. qs = a/273, cv = 1. (fi, = 1/12). 
4. A flat r2 m = 0.3 A-dominated universe. This flat universe with a cos- 
mological constant may be interpreted as a degenerate case of a solid dark 
matter component in which t s — > 1. 

Figures 1(a) and 1(b) indicate the evolution of the gravitational potentials <fi and 
if), respectively, in the limit k — > (i.e., on superhorizon scales) for the various models. 
The potentials have been normalized to unity at early times and the horizontal axis 
indicates conformal time, normalized so that rj — 1 today. Models 1, 2a-c, 3a-c, 4 are 
indicated, with no offset for model 1 and offsets increasing by .1 for each successive 
model, introduced to separate the curves in the plot. Although models (1) and models 
(2a)- (2c) have the same evolution of the scale factor, the evolution of the potential is 
different at later times. The fact that <fi and ip evolve differently is due to the presence 
of large anisotropic stresses. In the domain wall dominated models the decay of the 
Newtonian potential ip is much greater than in the hyperbolic or A models leading to 
a significantly larger integrated Sachs- Wolfe contribution to the CMB moments. 

Figure 2 shows the CMB moments for these models from i = 2 through £ = 30 
normalized so that C30 agrees for all the models. The vertical axis is q • £(£ + 1) 
with offsets increasing by of .2 The shapes of the CMB moments were computed 
with a method that does not take into account the effects producing the beginning 
of the rise toward the Doppler peak toward increasing I. In other words, l{£ + 1)q 
would be constant for a flat CDM universe when in fact there is a 40 % rise in this 
quantity by i = 30. Therefore, only the relative differences in shape are signficant. 
A more detailed study of these models using a Boltzman code will be presented in a 
forthcoming paper. 
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5. Discussion 



In this paper we have developed a continuum formalism for describing the dy- 
namics of a 'solid dark matter' (SDM) component and shown how cosmological per- 
turbations evolve with such a component included. The advantages of positing an 
SDM component are: (1). It is possible to reconcile a spatially flat cosmology with 
the numerous measurements of Q indicating a low value because most methods of 
measuring Q are sensitive only to matter that is clustered (e.g., on scales comparable 
to the size of galaxy clusters or smaller) and the SDM remains unclustered, except 
on the very largest scales comparable in size to our present horizon. (2). SDM can 
provide the negative pressure suggested by the recent SNIa observations at high red- 
shift, thus explaining the lower than expected apparent luminosities of the distant 
supernovae. Supernova observations can potentially constrain the equation of state, 
thus distinguishing a SDM-dominated universe from a cosmological constant domi- 
nated universe. (3). With SDM it is not necessary to posit a new, surprisingly small 
energy scale. SDM from string or domain wall networks results from new physics at 
higher energy scales. SDM thus avoids the fine tuning difficulties of a straightout A 
term or of the 'quintessence' models with an extremely slowly evolving scalar field 
that gives the same qualitative effect as A. 

Introducing an equation of state with negative pressure is a delicate matter. If 
one wishes to consider perturbations to homogeneity and isotropy, considerations of 
general covariance and causality prohibit one from introducing a smooth background 
component that does not cluster in an ad hoc way. Since locally it is impossible 
to determine what the 'unperturbed' background solution in the absence of pertur- 
bations would have been, a 'smooth component' must be introduced in a manner 
that specifies how perturbations evolve, and to do this more than merely specifying 
w = (p/p) is required. For an equation of state with negative pressure, to posit a 
perfect fluid is not allowed, because the sound speed on short wavelengths would be 
imaginary, indicating instabilities whose growth rate diverges as the wavelength ap- 
proaches zero. In SDM a sufficiently large shear modulus removes this instability. In 
the slowly rolling field models the instability is lacking for an entirely different reason: 
there is an inertia associated with changing the stress-energy. 

Physically, how the instability is avoided in SDM and in quintessence is manifested 
by the following qualitative differences: (1). SDM, unlike quintessence and most 
types of dark matter, generates anisotropic stresses. (2). SDM has vector modes 
with nonvanishing sound speed. (3). The resistance to pure shear in SDM gives the 
graviton a mass, changing the gravity wave contribution to the CMB on large angular 
scales. 
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We finally offer the following more technical remarks comparing SDM to other 
possible sources of negative pressure discussed in the literature: 



(1) . A pure cosmological constant may be regarded as a degenerate case of 
the action (2.3) with p{gm, <?(2), 9(3)) set to (A/Sn)y/g a i,. This degenerate case 
greatly enlarges the reparameterization invariance of (2.3). Because p = —p 
exactly, the smooth dark matter stress-energy no longer singles out a special 
time direction, and consequently for this special case general reparameteriza- 
tions that mix y and t are allowed. 

(2) . If p = 0, eqn. (2.3) becomes a Lagrangian description of a perfect fluid. 
In this special situation the Lagrangian description is much more cumbersome 
than the Eulerian description, especially with general relativity taken into ac- 
count. When p ^ 0, however, an Eulerian description is no longer possible. 
If p = and p = wp where w < and w / -1, the sound speed becomes 
imaginary, indicating an instability, most severe on the smallest scales. With- 
out resistance against pure shear, the solid dark matter component would 
be similarly unstable. However, for the non-Abelian strings and domain wall 
networks p is sufficiently large to stabilize the medium, as evidenced by all 
sound speeds being real. 

(3.) Quintessence differs in that the medium has an internal scalar degree of 
freedom. One could in fact write down a more general action for a low-energy 
description that combines quintessence with the continuous medium: 



S = 



dt dt (51) 



dt J d 3 y Vh p(g {l) ,g {2 ),9(3),4>)^-G. 
+ f d^XsJ\G^\ l - [(T? f 0) 2 - C 2(V, ^(1)^(2)^(3)) (2> y 0) 2 ] -V[<t>\. 

Because the stress-energy of the medium breaks Lorentz invariance down to 
the rotation group, c s need not equal the speed of light. 

4. Alexander Vilenkin' 291 has pointed out that for the special case w = 
— 1/3 (i.e., frustrated strings), due to a cancellation in Newtonian gauge the 
strings do not experience a gravitational force from nonrelativistic matter. In 
eqn. (2.21) this can be seen as a cancellation between the gradients of /ioo 
and ha that occurs only for this special case under the assumption that the 
source of these potentials has no anisotropic stresses. 
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Figure Captions 

Figure 1. Panels (a) and (b) show the evolution of the potentials and ip, respectively, 
for vanishing wavenumber as a function of conformal time for models 1, 2(a-c), 
3(a-c), and 4. All the potentials are normalized to unity at vanishing conformal 
time. The curves for each successive model have been displaced upward by 0.1. 

Figure 2. CMB Multipole Moments. The CMB multipole moments are plotted for models 
1, 2(a _ c), 3(a _ c), and 4 with successive curves displaced upward by 0.2. £(£+l)c£ 
is plotted and the curves are normalized to unity at £ = 30. 



24 



Figure 1 (a) 
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